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In the authors criticize our measurements [2] of the reduced two-point correlation 
function £(r) and of the one-point number variance cr 2 (r). On one hand, they confirm our 
conclusion that the estimator of £(r) for the initial conditions (IC) is in disagreement with 
the theoretical prediction at any scale in the simulation. However, the authors ascribe such a 
disagreement to the estimator's noise which they claim to be larger than the expected signal 
■ at any scale in the simulation. On the other hand, they find that the behavior of the estimator 
of <r 2 (r) is in good agreement with the theoretical expected behavior. Note that they do not 
provide any argument to explain the mismatch between the behaviours of £ (r) and a 2 (r) . In 
order to clarify the matter we have firstly recomputed both estimators and we have increased 
the statistical precision. Our result is that the new estimation of £(r) is in agreement with 
our previous one and with the one by pQ, while the new cr 2 (r) is in a better agreement with 
the theoretical behavior (see Fig.l). Hence, at least from the point of view of estimations, 
our conclusions are in good agreement with those of pQ. Our previous determination of o~ 2 (r) 
was affected by an imprecision in the normalization of the box which affected our estimation 
at scales comparable to the box size. The main point of disagreement with PP concerns the 
interpretation of the measured statistical estimators, and their relation with the expected 
theoretical behaviors. The following discussion will allow us to clarify also the mismatch in 
the determinations of cr 2 (r). 

The estimators of <J 2 (r) and £(r) are directly related (fluctuations in number of points in 
balls or shells centered on a random or distribution point normalized to the mean density [3]) 
and there is a precise theoretical relation with links these two statistics for a distribution with 
average density no: 

*»(r) = ^ \/x f v *vi(\*-m = + ^//* f v *vw-m, (i) 

where the second equality is explicitly satisfied by any discrete distribution; we have defined 
£(r) to be the non diagonal part of the complete correlation function £(r) = 5(r)/no + £(r) 
and hereafter we refer to the first term on the right hand side of Eq^ as the "discrete term" 
while the second part is denominated "correlated term" . Note that in the simulation we 
estimate the non diagonal part of the correlation function £s(r) (while the theoretical model 
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Fig. 1 - Normalized mass variance for the IC (Virgo) , theoretical behavior (CDM) and the contribution 
of the "discrete term" . The average distance between nearest neighbors (A) is also reported. 

gives the prediction for the complete £r(»") of a continuous field), and the estimator of <x 2 (r) 
contains both the discrete and correlated term (see |3] for more details). Now it is clear that, 
if the variance in the simulation cr|(r) agrees with the theoretical one cr^(r) we expect that 
£x(r) ~ £s( r ) if the discrete term in Eq^ is negligible with respect to the correlated term. 
However this is not the case: the contribution of 1/uqV is of order er^(r) at scales a few times 
larger than the mean interparticle distance (A) (see Fig.l). Hence our conclusion is that the 
agreement of the estimator of a 2 (r) with ct^ (r) is due to the fact that the particle distribution 
has still an average density no too small to have a mass variance dominated by the "correlated 
term" and discreteness effects are dominant at scale larger than l/10th of the box size. 

This is in agreement with the result that the behavior of the estimated £s (r) is different 
at all scales from the expected one of £t(t)' Its amplitude is very small on large scales and 
it becomes negligible with respect to discrete noise. A simple estimation of such a noise, by 
taking the lower limit of the Poisson case (i.e. the discrete term only), gives that <5£ ~ 1/ y/no ~ 
1/V256 3 « 2 x 10~ 4 (as V = 1 in our units). Using £t(j) one can see that |£r| < <5£ for r > 
0.05 w 10(A). Now, given the fact that at scales smaller than 3(A) w 0.015 £s is still oscillating 
between positive and negative values, due to the in-print of the pre-initial distribution, at best 
one can recover £t in the range w [3,10] x (A), while on larger scales the "discrete" noise 
predominates. Clearly such an agreement at intermediate scales is verified if there is not, in 
addition to "discrete" noise, a systematic difference between £x and £s: instead this seems to 
be the case as £s ^ £r at all scales. For the power spectrum (which we have not discussed 
in our paper) the situation is rather similar to the case of the mass variance as P(k)k 3 w 
<r 2 (r = 1/k) (but see We believe that a more complete theoretical understanding of 

the important problem of generating particle distribution with given correlations must be 
addressed by computing the effect of the application of a correlated displacement field on a 
generic particle distribution: We will discuss such a point in a forthcoming paper. 
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